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Resumen 


Las sequías meteorológicas son un fenómeno natural recurrente que origina una 
escasez de precipitación. La severidad de las sequías meteorológicas se estima 
a través de algoritmos establecidos, conocidos como índices de sequías. Uno de 
tales procedimientos, quizá el más simple, es el índice de reconocimiento de 
sequías O RDI (Reconnaissance Drought Index), que está basado en el cociente 
entre la precipitación y evapotranspiración potencial (ETP), ocurridas en un cierto 
lapso, seguido de meses. En este estudio se aplica el RDI en tres duraciones de 
sequía meteorológica, en cada una de las tres estaciones climatológicas 
seleccionadas de cada zona geográfica o climática del estado de San Luis Potosí, 
México, que fueron: Villa de Arriaga (Altiplano Potosino), Río Verde (Zona Media) 
y Xilitla (Región Huasteca). Los registros mensuales de precipitación y 
temperaturas media y mínima de cada estación abarcan más de 50 años. La ETP 
se estimó con cuatro métodos: 1) la fórmula de Penman-Monteith, que es el 
criterio de referencia; los criterios de 2) Thornthwaite; 3) Turc, y 4) Hargreaves- 
Samani. Los procedimientos operativos de estos criterios se exponen en los 
apéndices (ver más adelante). El análisis de los resultados indica que los RDI 
estimados con el método de Hargreaves-Samani es el que mejor reproduce los 
resultados de la fórmula de Penman-Monteith en las tres localidades climáticas 
procesadas. También el método de Turc conduce a resultados similares a los de 
referencia y por ello se puede establecer que el RDI es un índice de sequías 
robusto, que prácticamente no depende del método de estimación de la ETP. Al 
haber una diferencia notable en los procedimientos operativos de la fórmula de 
Penman-Monteith y del método de Hargreaves-Samani, este último es una 
solución práctica muy importante. 
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Palabras clave: sequías meteorológicas, evapotranspiración potencial, pruebas 
estadísticas, error cuadrático medio, error de sesgo medio, tipos de sequías 
meteorológicas (ligeras, moderadas, severas y extremas). 


Abstract 


Meteorological droughts are a recurring natural phenomenon that causes lack of 
precipitation. The severity of meteorological droughts is estimated by established 
algorithams known as drought indices. One such procedure, perhaps the simplest, 
is the Reconnaissance Drought Index, or RDI, which is based on the ratio 
between precipitation and potential evapotranspiration (PET) for a determined 
continuous period of months. In this study, the RDI is applied to three durations 
of meteorological drought at a weather station selected from each of the three 
geographic or climate zones in the state of San Luis Potosi, Mexico, which are: 
Villa de Arriaga (Potosino Plateau), Río Verde (Mean Zone), and Xilitla (Huasteca 
Region). The monthly rainfall records and average and minimum temperatures 
of each station cover more than 50 years. PET was estimated by four methods: 
(1) the Penman-Monteith formula, which is the reference method, (2) the 
Thornthwaite, (3) the Turc, and (4) the Hargreaves-Samani. The operating 
procedures for these criteria are detailed in appendices. The analysis of the 
results indicates that the RDIs estimated with the Hargreaves-Samani method 
are best for reproducing the results of the Penman-Monteith formula, in the three 
climatic locations processed. The Turc method also led to results similar to those 
of the reference. Therefore, it can be said that the RDI is a robust drought index, 
which practically does not depend on the method of estimating the PET. There is 
a noticeable difference in the operational procedures of the Penman-Monteith 
formula and the Hargreaves-Samani method. The latter is a practical solution 
that is worth mentioning. 


Keywords: Meteorological droughts, potential evapotranspiration, statistical 
tests, mean square error, mean bias error, types of meteorological drought 
(light, moderate, severe and extreme). 
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Introducción 


La sequía meteorológica es un fenómeno natural regional producido por la 
variabilidad del clima, el cual origina un decremento en la precipitación normal 
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de determinada área, durante un lapso importante y por ello tiene efectos 
adversos en la propia naturaleza y en la sociedad. La severidad de las sequías 
meteorológicas se estima por lo general con base en los llamados índices de 
sequía, los cuales varían en complejidad y van desde los que utilizan una sola 
variable climática como el SPI (Standardized Precipitation Index) hasta los que 
desarrollan un balance hídrico-edafológico como el PDSI (Palmer Drought 
Severity Index). Hao y Singh (2015) han encontrado que una sola variable no es 
suficiente para caracterizar las sequías debido a que son fenómenos naturales 
recurrentes originados por múltiples factores y por ello se han propuesto los 
índices de sequía multivariados. 


El enfoque de la variable latente permite el desarrollo de índices multivariados 
(Hao 8 Singh, 2015), consiste en establecer nuevas variables climáticas 
mediante una diferencia o cociente de otras variables de amplio significado físico, 
como la precipitación (P) y la evapotranspiración potencial (ETP) mensuales. 
Cuando se utilizó la diferencia P-ETP y el algoritmo operativo del SPI, se 
desarrolló el SPEI (Vicente-Serrano, Beguería, € López-Moreno, 2010), y cuando 
se empleó el cociente P/ETP se estableció el RDI (Reconnaissance Drought Index) 
o índice de reconocimiento de sequías, cuyo procedimiento operativo es bastante 
simple (Tsakiris 8 Vangelis, 2005; Tsakiris, Tigkas, Vangelis, 4 Pangalou, 2007; 
Vangelis, Tigkas, 8 Tsakiris, 2013). 


El objetivo de este estudio radica en exponer con detalle el proceso de cálculo 
del RDI anual, con duraciones de 3, 6 y 12 meses. Se aplica este procedimiento 
a los datos mensuales de precipitación, y de temperaturas medias y mínimas de 
tres localidades del estado de San Luis Potosí, representativas de sus tres zonas 
geográficas o climáticas. En este contraste del RDI, la ETP se estimó con cuatro 
métodos, cuya descripción detallada se expone en los apéndices, y que son: 1) 
la fórmula de Penman-Monteith, que fue el método de referencia; 2) el criterio 
de Thornthwaite; 3) el criterio de Turc, y 4) el criterio de Hargreaves-Samani. 
Se analizan los resultados y se formulan las conclusiones. 


Métodos y materiales 


Ecuaciones del RDI:.t 


El índice de reconocimiento de sequías (RDI) se calcula al inicio como el cociente 
entre la precipitación mensual acumulada y la respectiva evapotranspiración 
potencial, en los k meses considerados de cada año ¡ estudiado (Tsakiris 8 
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Vangelis, 2005; Tsakiris et al., 2007; Vangelis et al., 2013; Campos-Aranda, 
2014): 


k i 

2 j=1Pj 
k i 
Ea ET) 


(1) 


¡ASE 
Ak = 


En la ecuación anterior, k es la duración de la sequía meteorológica estudiada; 
j, el mes considerado, e ¡ varía de 1 a MA, que es el número de años del registro 
procesado (> 30). Como las magnitudes de al pueden ser representadas 
probabilísticamente por la distribución Log-normal; entonces los valores 
estandarizados del RDI se obtienen con facilidad mediante la siguiente ecuación: 


RDIy =%2 (2) 
y: 


en la cual: 
y = In(a,,) (3) 


En la ecuación (2), y es la media aritmética y oy, la desviación estándar de los 


valores yí. Los valores positivos del RDIst indican lapsos húmedos y los negativos 
son sequías meteorológicas, con la severidad siguiente: leves hasta -1.00; 
moderadas, fluctuando de -1.00 a -1.50; severas, variando de -1.50 a -2.00, y 
por último, extremas, menores que -2.00. Las duraciones comunes de k son 3, 
6, 9 y 12 meses; las tres primeras relativas a los meses seguidos con mayor 
porcentaje de precipitación, y la cuarta al lapso de enero a diciembre. Las 
duraciones menores a un año también pueden corresponder al periodo de 
crecimiento de los cultivos o a las épocas de demanda alta. Campos-Aranda 
(2014) expone un contraste entre el RDI:« y el SPI. 


Estimación de la ETP de referencia 


Hacia finales de la década de 1970, la Organización de las Naciones Unidas para 
la Agricultura y la Alimentación (FAO) formuló los lineamientos para la estimación 
de las demandas de agua de los cultivos (Doorenbos 4 Pruit, 1977). Avances en 
la investigación y evaluaciones más exactas sobre el uso del agua por los cultivos 
mostraron que el método de Penman, sugerido por la FAO, con frecuencia 
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sobreestima los requerimientos y que los criterios empíricos alternativos 
expuestos de manera variable se apegan a la realidad (Allen, Pereira, Raes, 8 
Smith, 1998). 


En mayo de 1990, la FAO organizó un panel de expertos e investigadores, en 
colaboración con la Comisión Internacional de Riego y Drenaje (ICID), y la 
Organización Meteorológica Mundial (WMO), para revisar los métodos de 
estimación de las demandas de los cultivos y establecer sus modificaciones. El 
panel recomendó la adopción de la fórmula de Penman-Monteith como método 
estándar para la estimación de la evapotranspiración de referencia (ETo) y 
aconsejó sobre los procedimientos de estimación de sus diversos parámetros 
meteorológicos (Allen et a/., 1998). 


En el Apéndice 1 se expone la fórmula de Penman-Monteith, incluyendo los 
procedimientos de estimación de sus parámetros basados en datos 
meteorológicos. En el Apéndice 2 se describen las recomendaciones de la FAO, 
para su aplicación con datos sólo climáticos. Esto último convierte a la fórmula 
de Penman-Monteith en un método aplicable y válido mundialmente para 
cálculos y contrastes de la ETo. Allen et al. (1998) indican que es preferible 
aplicar la fórmula de Penman-Monteith, incluso con las aproximaciones del 
Apéndice 2, que el uso de cualquier otro método empírico. 


Métodos empíricos de estimación de la ETP 


En el Apéndice 3, en sus ecuaciones (A.20) a (A.31), se exponen los 
procedimientos operativos de tres métodos empíricos de estimación de la 
evapotranspiración potencial mensual (ETP/), que son: el de Thornthwaite, el de 
Turc y el de Hargreaves-Samanli. 


Registros climáticos procesados 


Se puede considerar que el estado de San Luis Potosí está dividido en tres 
regiones climáticas: Altiplano Potosino, Zona Media y Región Huasteca. La 
primera tiene clima semiárido, la segunda templado-seco y la tercera cálido- 
húmedo. En cada una de estas regiones se buscaron las estaciones climatológicas 
de más amplio registro y con el menor número de datos mensuales faltantes de 
precipitación, y temperaturas media y mínima. Se seleccionaron las tres 
siguientes: Villa de Arriaga, Río Verde y Xilitla. En cada una de ellas, los escasos 
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datos faltantes de precipitación se consideraron igual a la moda, estimada con 
base en el ajuste de la distribución Gamma Mixta a todos los valores mensuales 
disponibles (Campos-Aranda, 2005a). Los datos faltantes de temperaturas media 
y mínima se estimaron con un procedimiento de interpolación que toma en 
cuenta la tendencia observada en tal mes, antes y después del valor faltante. 


Los valores promedio mensuales de precipitación (mm), y de temperaturas 
media y mínima (*C) de cada estación climatológica procesada se citan en la 
Tabla 1, así como sus respectivos periodos de registro. En la estación Villa de 
Arriaga, en el lapso de 2010 a 2014, se emplearon valores promedio mensuales 
de temperaturas debido a que sus registros disponibles abarcaron hasta el año 
2009. En la Figura 1 se muestra la ubicación de las tres estaciones climatológicas 
en el estado de San Luis Potosí. 


Tabla 1. Valores promedio mensuales de los elementos climáticos indicados en 
las tres estaciones climatológicas seleccionadas y procesadas del estado de San 
Luis Potosí, México. 


Des- Ene Feb Mar Abr Mayo | Junio | Julio | Ago Sep Oct Nov Dic Anual 
crip- 
ción 


Estación climatológica: Villa de Arriaga (longitud 101% 23" OG; latitud 21% 54' N; altitud 2 170 msnm; NA = 53 años) 


Precipi 13.0 7.6 7.0 10.4 30.9 57.2 72.1 56.8 62.7 25.4 5.6 8.9 357.7 
tación 
(1962- 
2014) 


Porcen 3.6 2.1 2.0 2.9 8.6 16.0 20.2 15.9 17.5 Lil 1.6 2.5 100.0 
taje de 
Precipi 
tación 


Temp. 13.0 13.9 15.9 19.2 20.9 20.8 19.8 19.5 18.8 16.6 14.6 13.4 17.2 
media 
(1962- 
2009) 


Temp. 4.3 4.9 7.0 9.5 112 | 114 | 113 | 11.0 | 10.5 8.4 6.0 4.8 8.4 
minim 
a 

(1962- 
2009) 


Radiac 400 380 505 650 650 550 520 420 380 450 430 390 480 
ión 
solar 


(cal/c 

m?/día 

) 

ETP: 84.1 87.8 126.5 160.6 | 177.7 157.1 151.6 133.6 114.3 115.5 97.8 88.7 1495. 
Penma 4 
n- 

Montei 

th 

ETP: 37.6 39.3 57.6 82.1 102.3 101.4 98.1 92.9 77.5 60.1 44.7 39.1 

Thornt 832.6 
hwaite 

ETP: 82.7 75.8 113.5 156.9 | 162.6 138.8 | 129.0 105.5 95.2 104.4 94.1 82.1 1340. 
Turc 7 
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ETP: 
Hargre 
aves- 
Saman 
¡ 


87.7 


77.4 


121.2 


166.5 


180.1 


146.9 | 139.9 


112.0 


96.1 


110.3 


96.0 


86.4 


1420. 


Esta 


ción clima 


tológica: 


Río Verde 


longitud 99% 59” OG; latitud 21% 56' N; 


altitud 987 msnm; NA = 54 años) 


Precipi 
tación 
(1961- 
2014) 


12.2 


10.8 


9.4 


32.7 


36.5 


88.7 


88.3 


71.7 


103.4 


44.2 


15.4 


12.9 


526.2 


Porcen 
taje de 
Precipi 
tación 


2.3 


2.0 


1.8 


6.2 


6.9 


16.9 


16.8 


13.6 


19.7 


8.4 


2.9 


2.5 


100.0 


Temp. 
media 
(1961- 
2014) 


16.2 


18.3 


21.7 


24.6 


26.4 


26.1 


25.0 


25.1 


23.9 


21.8 


19.0 


17.0 


22.1 


Temp. 
mínim 
a 
(1961- 
2014) 


8.6 


10.0 


12.6 


15.7 


18.2 


19.0 


18.3 


18.3 


17.7 


15.1 


12.1 


9.7 


14.6 


Radiac 
ión 
solar 


375 


350 


440 


550 


540 


540 


530 


480 


350 


390 


375 


350 


430 


Montei 
th 


78.4 


87.1 


130.6 


157.7 


170.0 


160.5 | 157.6 


148.1 


107.5 


110.3 


93.0 


84.7 


ETP: 
Thornt 
hwaite 


36.6 


47.5 


84.3 


119.3 


154.2 


147.0 | 135.4 


131.3 


107.1 


82.6 


54.3 


41.0 


ETP: 
Turc 


88.0 


81.2 


115.7 


149.0 


150.5 


149.8 | 145.0 


132.6 


98.2 


104.1 


95.0 


84.9 


ETP: 
Hargre 
aves- 
Saman 
¡ 


90.8 


124.4 


162.1 


171.9 


165.0 | 163.1 


147.8 


101.4 


110.5 


95.5 


Estación climatológica 


: Xilitla (longitud 98% 59' OG; 


latitud 21% 23'N; al 


titud 630 msnm; NA = 50 año 


Precipi 
tación 
(1965- 
2014) 


62.6 


65.3 


72.5 


115.3 


175.5 


373.9 | 432.2 


429.9 


566.1 


292.5 


101.5 


59.0 


Porcen 
taje de 
Precipi 
tación 


2.3 


2.4 


2.6 


4.2 


6.4 


13.6 


15.8 


15.7 


20.6 


10.7 


3.7 


2.1 


Temp. 
media 
(1965- 
2014) 


17.4 


18.7 


21.4 


24.2 


25.9 


26.2 


25.6 


25.9 


25.0 


23.1 


20.3 


18.3 


22.7 


Temp. 
mínim 


a 
(1965- 
2014) 


12.6 


13.4 


15.8 


18.5 


20.6 


20.9 


20.8 


20.2 


18.1 


15.6 


13.4 


17.6 


Radiac 
ión 
solar 
(cal/c 


350 


350 


400 


480 


500 


460 


490 


450 


310 


375 


360 


300 


400 
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m?/día 
) 
ETP: 69.1 78.3 106.1 127.9 | 142.9 132.1 | 137.9 | 132.2 98.4 105.9 87.4 74.3 1292. 
Penma 5 
n- 
Montei 
th 
ETP: 42.1 48.2 80.2 113.3 | 146.3 148.9 | 142.8 | 142.1 119.1 94.9 62.6 47.2 1187. 
Thornt 8 
hwaite 
ETP: 85.8 82.1 105.8 | 130.9 | 139.3 | 129.8 | 136.1 | 126.6 89.9 | 103.0 94.3 76.8 1300. 
Turc 3 
ETP: 87.9 82.5 112.4 | 140.2 | 157.2 141.1 | 152.9 | 141.3 92.2 110.0 95.1 77.2 1390. 
Hargre 0 
aves- 
Saman 
¡ 
1O2> 992 
240 —_—_—— 242 
Trópico de Cáncer 
San Luis Potosí ds: 
0 
22 222 
O Río Verde 
Villa de 
Arriaga Y 
Xilitla 
102% 99% 


Figura 1. Ubicación geográfica de las tres estaciones climatológicas procesadas 
del estado de San Luis Potosí, México. 


También en la Tabla 1 se exponen los valores promedio mensuales de la radiación 
solar (Rs;) obtenidos de los mapas de Hernández, Tejeda-Martínez y Reyes 
(1991), en sus páginas 65 a 76, para las tres localidades de las estaciones 
climatológicas seleccionadas y procesadas. Además, se muestran las magnitudes 
de la evapotranspiración potencial (mm) promedio mensual (ETP;), obtenidas a 
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través de las ecuaciones expuestas en los Apéndices 1 a 3, aplicadas a nivel 
mensual. 


Cuantificación de diferencias con el ECM y el ESM 


Con objeto de cuantificar las diferencias numéricas entre los valores anuales del 
RDIst de duración k debidas al impacto del cambio en el método de estimación 
de la ETPj se aplicaron los dos siguientes indicadores estadísticos: 1) el error 
cuadrático medio (ECM), y 2) error de sesgo medio (ESM), cuyas expresiones 
son (Vangelis et al., 2013): 


, a / 
ECM; = o NA (Xrefi - Xestk)'|] Ñ (4) 


1 . . 
ESM. => NA (Xrefi — Xesti) (5) 


En las expresiones anteriores, Xref; son los valores anuales del RDIst calculados 
con las ecuaciones (1) a (3), para una duración k en meses, con la ETP; estimada 
con la fórmula de Penman-Monteith (apéndices 1 y 2), que es el método de 
referencia, y Xestí son los mismos valores anuales del RDIst, pero calculados con 
base en una ETPj estimada con cada uno de los tres métodos empíricos 
expuestos (apéndice 3). Las ecuaciones (4) y (5) también se aplicaron a los 
valores anuales de la ETP' (k = 12). 


Análisis de los resultados 


Homogeneidad de los registros climáticos 


Con base en cada registro de precipitación, y temperatura media y mínima 
mensuales completado se obtuvieron sus valores anuales; con tales series se 
realizó el análisis de calidad estadística del registro, para lo cual se aplicaron 
siete pruebas: una general y seis específicas, las cuales fueron las siguientes: 1) 
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Von Neumann, que detecta pérdida de aleatoriedad por componentes 
determinísticas no especificadas; 2) Anderson, y 3) Sneyers, que buscan 
persistencia; (4) Kendall y 5) Spearman, que detectan tendencia; 6) Bartlett, 
que prueba la variabilidad, y 7) Cramer, que busca cambios en la media. En 
todas las pruebas se empleó un nivel de significancia (a) de 5%. Las pruebas 
estadísticas citadas se pueden consultar en WMO (1971), y Machiwal y Jha 
(2012). Los resultados de tales pruebas se han concentrado en la Tabla 2, donde 
NH y H significan serie o registro no homogéneo y homogéneo, respectivamente. 


Tabla 2. Resultados de la pruebas estadísticas aplicadas a los registros anuales 
de precipitación (P), y temperaturas media (Tt) y mínima (t) de las estaciones 
climatológicas indicadas. 


Villa de Arriaga Río Verde Xilitla 
Pruebas 
estadísticas P Tt t P Tt t P Tt t 
aplicadas 
1. Von Neumann NH NH NH NH NH NH H NH NH 
Zi Persistencia | NH NH NH NH NH NH H NH NH 
(Anderson) 
3. Persistencia | NH NH NH NH NH NH H NH NH 
(Sneyers) 
4, Tendencia (Kendall) H NH? H H NH? NH? H NH? 
5. Tendencia H NH? H H NH? | NH? H NH? 
(Spearman) 
6. Variabilidad H NH NH H H H H H H 
(Bartlett) 
7. Cambio en la media H NH H H NH NH H NH NH 
(Cramer) 
Coef. correlación lineal | 0.523 | 0.791 | 0.494 | 0.285 | 0.446 | 0.418 | 0.049 | 0.522 | 0.470 
de orden 1 (r:,) 


Con respecto a la precipitación anual (P), los registros de Villa de Arriaga y Río 
Verde muestran persistencia, detectada incluso con la prueba de Von Neumann. 
Al tomar en cuenta que la persistencia es una componente estadística de las 
series cronológicas, se puede continuar con los análisis encaminados a cuantificar 
sequías meteorológicas, pues los tres registros no presentan tendencia ni 
cambios en la media, es decir, pérdida de homogeneidad. 


Lo contrario ocurre con los registros de temperatura media anual (7t), los cuales 
son no homogéneos, ya que presentan persistencia (asociada con r,), tendencia 
ascendente y cambio en la media. En relación con los registros de temperatura 
mínima anual (t), el de la estación Río Verde es el menos homogéneo, pues 
presenta tendencia ascendente, y el de Villa de Arriaga es el más homogéneo, 
ya que sólo tiene persistencia. El de Xilitla también es no homogéneo, al tener 
persistencia y cambio en la media. 
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Resultados numéricos del RDIst 


Con base en los porcentajes de la precipitación mensual expuestos en la Tabla 
1, se definió el lapso de tres meses con más lluvia de julio a septiembre, y el de 
seis meses de mayo a octubre. En las estaciones de Villa de Arriaga, Río Verde 
y Xilitla, las sumas de porcentajes para k = 3 meses fluctúan en el 52%, y para 
k = 6 meses en 83%. Definidos tales lapsos se aplicaron las ecuaciones 1 a 3 a 
los datos de precipitación y a las estimaciones de la ETP; únicamente parte de 
tales resultados se muestran en las Tablas 3 a 5 debido a limitaciones de espacio. 


Tabla 3. Datos, valores del RDI:: y tipos de sequías meteorológicas anuales 
(TS) calculados con la evapotranspiración potencial estimada según los criterios 
empíricos indicados en la estación climatológica Villa de Arriaga, San Luis 


Potosí. 
p Según fórmula de Penman-Monteith Según método de Hargreaves-Samani 
Año ETP” | k= 3 meses k=6 k=12 ETP k=3 k=6 k=12 
(mm) meses meses meses meses meses 
(mm) (mm) 

RDIst TS RDIs | TS RDIst | TS RDIs | TS RDIs | TS RDls | TS 

1962 268.4 1495.8 | -0.110 SsL a SL E SL | 1433.6 + SL iS SsL a SsL 
0.169 0.217 0.112 0.179 0.229 

1963 279.1 1442.0 | -0.444 SsL de SsL a SL | 1387.2 s SL a SsL S SsL 
0.755 0.089 0.474 0.786 0.108 

1964 489.0 1392.8 0.600 el 2 $ 1370.2 = > ES 
0.960 0.917 0.553 0.889 0.855 

1965 426.1 1376.3 0.600 a 2 a 1364.0 E > > 
0.389 0.704 0.551 0.306 0.632 

1966 631.7 1293.4 el > E 1287.1 E > ES 
0.5714 1.264 1.475 0.491 1.165 1.391 

1967 538.0 1299.7 0.722 cs > a 1275.1 de > > 
0.991 1.195 0.668 0.913 1.137 

1968 379.9 1368.5 0.664 = = > 1347.8 ie - S= 
0.584 0.520 0.595 0.500 0.459 

1969 174.0 1450.0 | -0.285 SsL 3 SsL E SL | 1405.9 iS SsL s SsL - SsL 
0.786 0.897 0.331 0.841 0.924 

1970 269.7 1418.5 | -0.078 SL se a SL | 1381.9 > SsL iS SsL a SL 
0.045 0.119 0.133 0.009 0.159 

1971 593.0 1394.1 1.259 se a e 1372.1 - = se 
1.577 1.241 1.191 1.489 1.177 

1972 347.0 1417.3 0.248 > SS > 1388.5 Ñ S se 
0.603 0.308 0.194 0.533 0.257 

1973 328.5 1418.9 0.647 $e 2 e 1390.9 - - - 
0.538 0.213 0.586 0.458 0.162 

1974 156.9 1423.4 0.021 pe = SsL - SM | 1391.7 ie SsL - SM - SM 
0.978 1.041 0.031 1.044 1.081 
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1975 280.5 | 1432.3 0.329 S a 3 SL | 1397.9 a - e SL 
0.185 0.069 0.267 0.118 0.112 
1976 359.5 | 1397.9 0.695 y a S 1362.3 - - > 
0.741 0.391 0.638 0.669 0.348 
1977 96.0 1395.6 | -0.319 | SL - ss - SS | 1353.3 E SL 5 ss E ss 
1.535 1.837 0.338 1.567 1.860 
1978 507.5 | 1363.3 1.321 S - S 1325.2 - - > 
1.333 1.016 1.288 1.273 0.974 
1979 192.0 | 1419.7 0.077 S - SL y SL | 1372.7 y a SL 5 SL 
0.770 0.695 0.051 0.805 0.719 
1980 375.5 | 1393.4 0.064 y - SL S 1355.5 - E SL a 
0.062 0.470 0.039 0.101 0.430 
1981 270.9 | 1385.6 | -0.567 | SL E SL = SL | 1340.1 = SL = SL - SL 
0.343 0.072 0.577 0.367 0.100 
1982 292.5 | 1438.5 | -0.443 | SL y SL a SL | 1376.0 = SL = SL = SL 
0.015 0.006 0.449 0.024 0.015 
1983 154.5 | 1385.6 | -0.219 | SL E SL E SM | 1345.3 = SL = SL = SM 
0.791 1.021 0.232 0.804 1.050 
1984 148.0 | 1419.9 | -0.017 | SL E SM ñ SM | 1372.2 = SL = SM = SM 
1.055 1.135 0.059 1.090 1.156 
1985 167.0 | 1428.0 | -0.765 | SL - SL a SL | 1369.2 = SL = SL = SL 
0.792 0.941 0.760 0.789 0.949 
1986 504.0 | 1404.8 | -0.621 | SL a S 1371.9 E SL - E 
1.012 0.953 0.638 0.978 0.904 
1987 504.9 | 1377.7 0.939 S a S 1321.5 S - = 
1.049 0.989 0.947 1.056 0.970 
1988 256.0 | 1407.2 0.337 = S ñ SL | 1343.9 s 0.51 - 5 SL 
0.054 0.194 0.338 0.200 
1989 359.5 | 1469.2 | -0.052 | SL - SL - 1354.1 E - SL - 
0.155 0.307 0.018 0.074 0.358 
1990 | 728.3 | 1507.6 0.916 y a - 1384.9 - - - 
1.440 1.456 0.989 1.518 1.507 
1991 | 1028.5 | 1564.4 1.572 a y - 1416.6 y S = 
2.045 1.977 1.649 2.138 2.049 
1992 | 838.0 | 1624.5 0.015 a - - 1458.3 - S = 
0.934 1.567 0.115 1.051 1.656 
1993 523.0 | 1716.8 0.696 : 3 - 1639.7 S - - 
0.754 0.677 0.634 0.690 0.667 
1994 564.9 | 1561.4 0.632 y 3 - 1448.4 - - - 
1.123 0.968 0.706 1.191 1.005 
1995 392.8 | 1530.0 | -0.103 | SL - SL - 1409.2 a SL - SL - 
0.190 0.388 0.042 0.115 0.440 
1996 | 666.0 | 1686.9 1.263 7 - = 1506.8 - - E 
1.136 1.115 1.370 1.260 1.215 
1997 | 454.0 | 1730.8 | -0.052 | SL 3 - 1536.6 Ss - - 
0.272 0.424 0.043 0.393 0.538 
1998 | 418.0 | 1748.9 0.181 3 a - 1564.4 S - - 
0.383 0.267 0.285 0.504 0.369 
1999 249.0 | 1698.0 | -0.126 | SL = SL = SL | 1516.1 = SL = SL = SL 
0.256 0.558 0.065 0.153 0.449 
2000 123.0 | 1599.0 | -5.157 | SE = SM + SS | 1469.4 a SE 5 SM E ss 
1.464 1.648 5.148 1.411 1.582 


158 


2001 111.8 | 1564.2 | -0.470 | SL - SM - SS | 1458.5 E SL a SM = ss 
1.484 1.773 0.452 1.438 1.730 
2002 139.3 | 1582.7 | -0.956 | SL a SM - SM | 1516.8 5 SL = SM a SM 
1.437 1.421 0.956 1.419 1.426 
2003 337.0 | 1638.4 0.463 S y = 1590.3 E - E SL 
0.356 0.014 0.455 0.343 0.021 
2004 | 613.4 | 1645.2 0.452 :S a S 1532.7 - - o 
0.978 1.019 0.481 1.017 1.048 
2005 287.0 | 1667.3 0.147 = - SL - SL | 1548.7 S a SL S SL 
0.323 0.287 0.195 0.262 0.246 
2006 340.0 | 1637.5 0.120 : a S 1531.1 - - S 
0.121 0.029 0.152 0.158 0.058 
2007 270.0 | 1630.7 | -0.564 | SL E SL A SL | 1531.0 S SL = SL = SL 
0.473 0.353 0.537 0.439 0.329 
2008 546.4 | 1602.8 1.214 : S = 1499.5 - - a 
1.177 0.867 1.247 1.207 0.890 
2009 201.0 | 1518.3 | -0.501 | SL > SL E SL | 1468.8 = SL = SL SL 
0.627 0.731 0.497 0.618 0.755 
2010 264.0 | 1498.3 0.029 S E SL E SL | 1420.1 - = SL = SL 
0.763 0.248 0.036 0.753 0.241 
2011 75.8 1498.3 | -1.324 | SM - SE ñ SE | 1420.1 = SM = SE = SE 
2.076 2.357 1.323 2.065 2.338 
2012 134.9 | 1498.3 | -0.669 | SL - ss ña SM | 1420.1 = SL = ss = SM 
1.610 1.382 0.665 1.599 1.369 
2013 138.9 | 1498.3 | -1.436 | SM - ss a SM | 1420.1 - SM = ss = SM 
1.849 1.333 1.436 1.838 1.320 
2014 165.0 | 1498.3 | -1.518 | SS A SM E SM | 1420.1 = ss = SM = SM 
1.287 1.043 1.517 1.276 1.030 


Tabla 4. Datos, valores del RDIs: y tipos de sequías meteorológicas anuales 
(TS) calculados con la evapotranspiración potencial estimada según los criterios 
empíricos indicados en la estación climatológica Villa de Arriaga, San Luis 


Potosí. 
Te e Según criterio de Thornthwaite Según criterio de Turc 
o esoo PEBTE k=3 k=6 k= 12 ETP k=3 k=6 k= 12 
( ( meses meses meses meses meses meses 
) ) (mm) (mm) 
RDIst | TS | RDÍst | TS | RDIst | TS RDIst | TS | RDÍst | TS | RDIst | TS 
196 17. 9.0 5 SL a sL a SL | 1354.2 > SsL = SL + SsL 
2 4 825.3 0.12 0.20 0.21 0.11 0.17 0.22 
2 3 5 2 5 8 
196 16. 8.3 776. E SsL £ SsL E SL | 1319.6 + SL + SsL < SsL 
3 5 0.38 0.70 0.04 0.48 0.79 0.12 
3 2 5 5 9 0 
196 16. 9.0 - - - 1309.5 a = Ñ 
4 2 758.4 0.64 1.10 0.94 0.53 0.85 0.83 
4 7 6 8 9 0 
196 15. 8.9 - - - 1298.7 Al Al a 
5 9 756.8 0.62 0.41 0.71 0.53 0.28 0.61 
3 6 6 7 5 4 
196 14. 7.5 A E Ñ 1217.5 Ñ - Ñ 
6 0 700.3 0.63 1.33 1.51 0.47 1.13 1.37 
0 5 5 6 7 9 


196 | 13. 6.6 - o 3 1206.2 - y - 
7 7 693.0 | 0.76 1.08 1.26 0.65 0.88 1.12 
6 9 0 1 6 6 
196 | 15. 8.3 E y e 1282.4 - - y 
8 5 742.8 | 0.70 0.63 0.55 0.57 0.48 0.44 
9 9 3 9 0 3 
196 | 16. 8.8 y SL | 0.85 | SL | 0.89 | SL | 1330.8 y SL - SL - SL 
9 8 798.5 | 0.30 2 4 0.33 0.84 0.92 
1 4 0 3 
197 | 16. 8.5 s SL E = SL | 1312.1 S SL S SL S SL 
0 2 772.4 | 0.08 0.07 0.09 0.14 0.02 0.16 
6 d: 4 3 6 7 
197 | 16. 8.8 y a E 1309.9 y - y 
1 1 760.1 | 1.31 1.65 1.26 1.17 1.46 1.15 
0 8 9 2 0 1 
197 | 16. 8.9 - - - 1321.0 S A S 
2 4 774.1 | 0.28 0.66 0.32 0.18 0.51 0.24 
4 4 9 1 1 2 
197 | 16. 9.0 - - - 1323.0 y S y 
3 5 778.0 | 0.66 0.55 0.22 0.57 0.43 0.14 
9 5 8 1 8 8 
197 | 16. 9.0 - - SL - Ss 1325.1 y SL a Ss 1.08 Ss 
4 6 780.5 | 0.00 0.98 1.03 M 0.03 1.05 M 8 M 
2 1 1 9 0 
197 | 16. 9.0 - y a SL | 1331.4 - - E SL 
5 6 780.2 | 0.42 0.27 0.04 0.25 0.10 0.12 
3 1 5 1 1 6 
197 | 15. 8.3 - - - 1300.8 - - - 
6 9 751.8 | 0.76 0.84 0.43 0.62 0.64 0.32 
5 8 9 0 2 7 
197 | 15. 7.9 y SL - Ss 1.79 | SS | 1288.2 : SL = SS | 1.86 | SS 
7: 6 747.4 | 0.25 1.47 M 0 0.35 1.57 1 
7 5 1 8 
197 | 14. 7.3 y e S 1259.3 S - E 
8 9 727.8 | 1.36 1.42 1.07 1.27 1.24 0.95 
> 8 8 0 4 7 
197 | 16. 8.2 0.92 | - E SL = SL | 1301.9 > S SL F SL 
9 0 768.0 0.74 0.66 0.04 0.81 0.72 
1 1 4 4 2 
198 | 15. 8.2 S E SL S 1287.1 0.29 | - y SL > 
0 7 754.5 | 0.05 0.05 0.50 0.11 0.41 
3 4 7 7 7 
198 | 15. 7.6 5 SL 5 SL = SL | 1276.0 = SL 7 SL + SL 
1 4 737.9 | 0.50 0.25 0.00 0.58 0.38 0.11 
3 3 9 9 8 3 
198 | 16. 7.9 y SL - - 1308.4 y SL a SL o SL 
2 1 766.2 | 0.37 0.06 0.05 0.46 0.04 0.02 
1 7 7 0 1 7 
198 | 15. 7.7 - SL A SL a SL | 1279.0 E SL E SL a Ss 
3 4 739.9 | 0.08 0.67 0.96 0.24 0.81 1.05 M 
9 6 6 8 6 5 
198 | 16. 8.2 S a SL A Ss 1307.6 a SL - >: z >: 
4 1 762.0 | 0.04 0.96 1.08 M 0.07 1.10 M 1.16 M 
3 8 9 3 3 4 
198 | 16. 7.8 - SL a SL 5 SL | 1301.6 - SL - SL - SL 
5 0 764.1 | 0.72 0.75 0.88 0.76 0.79 0.95 
4 4 9 6 8 4 


198 | 16. | 8.8 - [sL a - [| 1303.5 2 SE E S 
6 2 767.6 | 0.62 1.05 0.97 0.64 0.95 0.88 
7 2 7 4 9 8 
198 | 14. | 6.8 E . - | 12447 - E 
7 8 732.4 | 0.94 1.08 1.05 0.93 1.03 0.96 
4 4 9 3 4 8 
198 | 15. | 6.9 ] . - [sL| 1269.2 - - SE 
8 2 745.8 | 0.36 0.08 0.12 0.32 0.03 0.19 
4 8 3 7 9 9 
198 | 15. | 5.9 - Isi| - |sL - | 1281.7 - - [SL A 
9 6 755.8 | 0.01 0.05 0.43 0.01 0.08 0.35 
4 6 0 2 4 2 
199 | 16. | 6.1 E E - | 12998 ] a E 
0 e, 790.6 | 0.92 1.47 1.55 0.98 1.50 1.50 
6 2 1 5 9 7 
199 | 16. | 5.7 - s - | 13275 - - - 
1 9 821.0 | 1.62 2.08 2.07 1.64 2.13 2.04 
1 5 2 5 7 8 
199 | 17. | 6.2 E - - | 1362.0 - - - 
2 9 881.9 | 0.09 0.85 1.60 0.14 1.07 1.66 
2 5 3 2 E 3 
199 | 22. | 12. | 1603. -= [sii - |sL|o0.21|sL| 1454.2 - - - 
3 2 8 7 0.37 0.58 0 0.78 0.84 0.76 
9 8 6 5 7 
199 | 17. | 7.6 - E - | 13457 E E E 
4 6 872.3 | 0.56 0.97 0.95 2 0.72 1.21 1:02 
4 7 3 2 6 5 
199 | 16. | 6.3 2 (Sl st - | 1306.3 ss sE E 
5 5 821.4 | 0.00 0.26 0.43 0.04 0.07 0.46 
2 0 9 4 9 8 
199 | 19. | 6.9 E E - | 1387.1 - - - 
6 2 981.5 | 1.12 0.83 1.03 1.40 it 1.24 
6 7 2 ] 5 9 
199 7.0 - [sL| 041] - - | 1429.4 - - - 
O E 996.7 | 0.24 0.35 0.87 0.44 0.55 
9 2 6 4 3 9 
199 | 20. | 8.3 | 1055. sE E - | 1451.1 - a - 
8 7 0 0.11 0.05 0.12 0.34 0.56 0.39 
2 6 0 1 1 6 
199 | 19. | 6.7 - [sii - |sL 0.53 | sL | 1420.8 = [sti - [sil - |sL 
9 3 923.3 | 0.06 0.23 2 0.06 0.13 0.43 
0 4 1 2 4 
200 | 18. | 7.7 - [sel - Ss - [| ss| 1389.3 - [seE| - 3 - [ss 
0 3 854.6 | 5.14 1.38 | M | 1.59 5.13 1.40 | M | 1.57 
4 5 7 3 2 4 
200 | 18. | 8.2 A Ss - [ss| 1380.1 A E E Ss - [ss 
1 1 845.1 | 0.40 1.44 | M | 1.74 0.45 1.43 | M | 1.72 
7 6 0 6 0 2 
200 | 19. | 10. a S . Ss | 1429.5 - [sil - Ss E 3 
2 6 7 | 912.0 | 0.92 1.42 | M | 149 | M 0.95 1.40 | M| 141 | M 
0 2 6 2 2 4 
200 | 21. | 13. | 1019. a S - [| sSL| 1486.5 E E - [| sL 
3 5 1 0 0.48 0.34 0.18 0.46 0.35 0.00 
8 6 6 1 9 4 
200 | 20. | 10. a a - | 1446.7 E E E 
4 2 1 | 927.9 | 0.55 1.14 0.98 0.48 1.01 1.04 
5 8 8 1 sb 2 


200 | 20. 10. - - SL - SL | 1458.5 - a SL - SL 
5 6 3 951.2 | 0.23 0.18 0.34 0.19 0.25 0.24 
4 2 2 8 8 0 
200 | 20. 10. E y S SL | 1446.0 - - y 
6 2 3 924.2 | 0.22 0.28 0.00 0.15 0.15 0.05 
8 6 6 1 5 7 
200 | 20. 10. in SL E SL dE SL | 1446.0 S SL - SL - SL 
7 2 5 925.3 | 0.48 0.35 0.39 0.53 0.43 0.32 
5 0 9 4 6 8 
200 | 19. 9.8 y E E 1409.1 y y e 
8 3 901.2 | 1.28 1.27 0.84 1.24 1.19 0.89 
9 8 sl 4 9 E 
200 | 18. 10. a SL e SL - SL | 1386.1 y SL a SL - SL 
9 4 2 861.4 | 0.49 0.60 0.77 0.49 0.61 0.75 
8 6 8 6 6 0 
201 17. 8.4 - - SL + SL | 1348.4 a = SL - SL 
0 2 808.6 | 0.04 0.74 0.20 0.03 0.75 0.24 
7 5 8 1 5 9 
201 17. 8.4 - Ss - SE - SE | 1348.4 a Ss E SE a SE 
Bl 2 808.6 | 1.32 M 2.08 2.32 1.32 M 2.05 2.33 
6 3 4 El 7 2 
201 17. 8.4 - SL - SS E Ss 1348.4 3 SL e ss Ss Ss 
2 2 808.6 | 0.66 1.60 1.32 M 0.66 1.59 1.37 M 
1 8 4 8 5 0 
201 17. 8.4 - Ss - SE y Ss 1348.4 S S: Sl ss y S: 
3 2 808.6 | 1.44 M 1.85 1.29 M 1.43 M 1.83 1.32 M 
0 2 7 5 2 1 
201 17. 8.4 - ss - Ss - Ss 1348.4 - ss - Ss - Ss 
4 2 808.6 | 1.52 1.27 M 1.00 M 1.51 1.27 M 1.03 M 
3 9 5 7 4 3 


Tabla 5. Datos, valores del RDI:: y tipos de sequías meteorológicas anuales 
(TS) calculados con la evapotranspiración potencial estimada según los criterios 
empíricos indicados en la estación climatológica Xilitla, San Luis Potosí. 


p Según fórmula de Penman-Monteith Según método de Turc 
Año ETP' k=3 k=6 k = 12 ETP k=3 k=6 k = 12 
(mi) meses meses meses meses meses meses 
(mm) (mm) 
RDIst | TS | RDIs | TS | RDlst | TS RDIst | TS RDIst TS | RDIst | TS 
196 3010. 1248. - - a 1293. Ñ 0.517 - - 
5 6 6 0.75 0.61 0.60 6 0.69 0.52 
5 6 6 2 2 
196 2627 1228. - SE - sL Ñ 1272. Ñ SE | -0.214 | SL - 
6 7 8 2.20 0.10 0.12 7 2.31 0.00 
9 6 3 1 9 
196 2912 1280. - - Sl 1290. - 0.501 - - 
7 1 6 0.64 0.44 0.37 1 0.60 0.39 
0 6 0 0 1 
196 2302. 1268. 2 - sL - sL 1289. = -0.345 | SL 5 SL 
8 7 2 0.41 0.28 0.53 1 0.34 0.61 
0 8 6 8 2 
196 3532. 1285. - - - 1301. - 1.172 a a 
9 0 4 1.44 1.10 1.13 0 1.50 1.18 
8 2 3 4 2 
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197 | 2961. | 1244. - - = "1273. - | 0.758 | - - 
0 0 1 0.22 0.77 0.55 4 0.18 0.51 

7 9 4 9 8 
197 | 3330. | 1301. z E - | 1300. - | 0.951 | - - 
1 5 2 0.74 0.94 0.84 9 0.73 0.93 

9 5 7 0 1 
197 | 3202. | 1238. - - - | 1292. - [SL | 0.639 | - - 
2 0 5 0.07 0.77 0.88 6 0.01 0.78 

2 1 8 1 9 
197 | 3166. | 1215 , S - | 1290. Sos | - 
3 0 6 0.27 1.05 0.91 1 0.12 0.74 

6 4 7 1 9 
197 | 3044. | 1247. - - - | 1290. - | 0.579 | - - 
4 9 6 0.86 0.71 0.65 9 0.75 0.57 

5 0 5 4 9 
197 | 3474. | 1294. - - - | 1294. - | 1.377 | - - 
5 5 7 1.30 1.18 1.03 0 1.51 1.13 

4 3 8 0 5 
197 | 3297. | 1167. - - - | 1263. - [| 0.828 | - - 
6 2 0 1.07 1.01 1.24 á 0.94 1.01 

3 4 5 8 4 
197 | 1787. | 1255. -= |ss| - |ss| - [ss] 1285. - (SE|-1656|ss| - |ss 
7 2 0 1.91 1.60 1.51 5 2.00 1.68 

3 1 6 8 6 
197 | 3276. | 1234. E - - | 1287. - | 0.886 | - - 
8 6 0 0.78 0.90 0.99 7 0.76 0.90 

1 9 5 1 4 
197 | 2364. | 1326. -= si] - [|sLÍ| - sL | 1297. - [sL|-0622|sL| - |sL 
9 1 2 0.11 0.69 0.61 6 0.02 0.52 

1 2 1 5 7 
198 | 1896. | 1363. = [sip - [ss| - [ss] 1314. - [SL - ss| - |ss 
0 6 4 0.61 1.79 1.61 2 0.44 14.70 1.52 

7 0 0 4 4 7 
198 | 3351. | 1307. - - - | 1304. - | 0.713 | - - 
1 7 4 0.18 0.62 0.85 6 0.24 0.94 

3 El 3 9 6 
198 | 1951. | 1345. -= [ss - [|ssÍ| - s | 1310. - (SS |-1.752|SS| - Ss 
2 0 8 1.81 1.77 1.44 | M 0 1.66 1.39 | M 

2 0 4 7 > 
198 | 3728. | 1320. - - - | 1296. - | 1.561 | - - 
3 5 8 1.79 1.39 1.24 Z 1.90 1.42 

0 1 2 9 8 
198 | 3758. | 1257. - - - | 1287. - | 1.718 | - - 
4 7 0 1.47 1.72 1.47 5 1.40 1.49 

3 0 4 0 3 
198 | 2720. | 1296. -= si - [| sL - | 1294. - [| sL|-0.222 | SL - 
5 6 9 0.44 0.17 0.04 7 0.48 0.08 

9 9 5 2 4 
198 | 2552. | 1284. - s - [si | - sL | 1299. - s |-0.299 | sL| - ]|sL 
6 2 9 1.433 | M | 0.24 0.17 3 1.44 | M 0.20 

3 4 5 9 5 
198 | 2388. | 1257. - Poo7 [si | - SL. (1272, - | -0.097 | sL| - | sL 
7 8 9 0.31 6 0.35 9 0.31 0.40 

9 6 0 1 
198 | 2815. | 1253. - | sL - - | 1288. - [| sL| 0.123 | - - 
8 8 0 0.16 0.22 0.32 9 0.19 0.25 

7 7 3 8 1 


198 2713 1333. S SL = SL 2 SL | 1303. - SL | -0.859 | SL a 
9 6 0 0.38 0.85 0.07 6 0.42 0.04 
6 0 6 3 4 
199 2639. 1285. - = SL = SL | 1302. E -0.291 | SL > SL 
0 0 6 0.29 0.23 0.04 9 0.23 0.07 
0 1 2 2 4 
199 3597. 1256. - - - 1301. - 1.194 - - 
1 0 3 0.85 1.19 1.29 9 0.77 1.25 
3 6 9 9 7 
199 3175. 1191. E y e 1279. - 0.451 - - 
2 D: 0 0.32 0.58 1.01 9 0.23 0.79 
0 5 2 6 6 
199 3481 1240. - - a 1296. S 0.922 - - 
3 9 3 0.55 1.01 1.21 0 0.47 1.13 
7 0 9 7 7 
199 2532. 1271. - - - SL | 1305. - -0.019 | SL 3 SL 
4 3 3 0.33 0.02 0.16 9 0.35 0.26 
3 9 4 1 0 
199 2617 1292. - - - SL | 1313 - 0.009 - - SL 
5 0 1 0.36 0.00 0.09 3 0.35 0.14 
2 2 6 6 3 
199 1918. 1305. - SL - Ss - Ss 1304. - SL | -1.107 | S - Ss 
6 0 1 0.69 1.08 M 1.38 M 4 0.71 M 1.44 M 
5 6 9 3 6 
199 1999. 1273. - SE - Ss - Ss 1302. - SE | -1.344 | S - Ss 
7 5 9 2.12 1.26 M 1.12 M 9 2.15 M 1.26 M 
6 9 3 6 3 
199 2819 1314. - - - 1326. - 0.201 y - 
8 5 2 0.56 0.15 0.13 8 0.57 0.13 
7 5 6 2 2 
199 2576. 1312. E a = sL | 1313 S 0.198 - = SL 
9 1 0 0.72 0.18 0.22 4 0.71 0.21 
9 8 1 0 1 
200 1974. 1307. E ss y Ss - Ss 1317. a SS | -1.309 | S - Ss 
0 6 4 1.86 1.21 M 1.27 M 6 1.89 M 1.36 M 
0 4 9 8 5 
200 2751 1293. e SL - SL e 1313 - SL | -0.059 | SL - 
1 0 4 0.15 0.05 0.10 9 0.16 0.06 
8 3 1 1 9 
200 1880. 1314. - SL - Ss - Ss 1309. - SL | -0.996 | SL - ss 
2 7 4 0.90 1.04 M 1.49 M 7 0.88 1.54 
2 0 7 8 8 
200 2686. 1294. - - - 1305. - 0.187 - E SL 
3 6 2 0.55 0.20 0.00 0 0.44 0.00 
Z 4 3 8 4 
200 2155. 1283. : SL a SL 7 SL | 1301. - SL | -1.041 Ss - SL 
4 8 1 0.75 0.96 0.84 5 0.76 M 0.93 
4 7 9 1 6 
200 2517 1332. - SL y SL - SL | 1316. - SL | -0.076 | SL = SL 
5 9 2 0.60 0.15 0.37 0 0.50 0.31 
1 2 5 7 8 
200 1654. 1379. - Ss - SE a SE | 1322. - S | -2.422 | SE - SE 
6 4 0 1.44 M 2.44 2.20 8 1.37 M 2.14 
3 8 6 8 0 
200 2948. 1311 a S S 1305. - 0.499 y a 
7 9 1 0.42 0.41 0.32 6 0.48 0.39 
0 0 6 6 4 


200 3475. 1332. S E - 1310. - 1.247 y a 
8 2 7 0.94 1.14 0.92 1 0.94 1.08 
3 3 2 0 3 
200 2239. 1393. a SL > SL 5 Ss 1325. - SL | -0.688 | SL = SL 
9 0 9 0.89 0.88 1.03 M 3 0.71 0.85 
2 7 0 4 1 
201 2880. 1340 - - SL e 1302. - 0.102 - - 
0 7 4 0.71 0.01 0.14 2 0.80 0.30 
0 6 3 5 4 
201 1554 1387 - Ss - SE - SE | 1313. - S | -2.045 | SE = SE 
1 5 8 1.16 M 2.12 2.48 4 1.10 M 2.37 
9 2 3 7 6 
201 2202. 1376. - SL - SL - Ss 1315. s SL | -0.905 | SL - SL 
2 8 5 0.19 0.98 1.04 M 2 0.10 0.88 
3 7 6 8 8 
201 3764. 1327. - - - 1304. A 1.297 - - 
3 8 5 1.33 1.18 1.26 0 1.41 1.44 
3 6 0 0 6 
201 3100. 1356 - SL - - 1311. - SL | 0.530 - - 
4 4 1 0.45 0.46 0.39 7 0.41 0.58 
2 4 2 8 8 


En la Tabla 3 y Tabla 4 se exponen todos los resultados de la estación Villa de 
Arriaga, la cual presenta mayor variabilidad, pues su precipitación y temperatura 
media anuales fluctúan de 75.8 a 1 028.5 mm y de 13.7 a 21.5 *“C. En la Tabla 
5 se expone una parte de los resultados correspondientes a la estación Xilitla, la 
de menor dispersión, como se observa en sus siguientes intervalos de fluctuación 
de precipitación y temperatura media anuales: de 1 554.5 a 3 764.8 mm y de 
21:2a:239 90€, 


En la Tabla 3, Tabla 4 y Tabla 5 se emplea la simbología siguiente para la 
severidad o tipos de sequía meteorológica: sequías ligeras (SL), sequías 
moderadas (SM), sequías severas (SS) y sequías extremas (SE). 


Los resultados numéricos del RDIst mostrados en las Tablas 3 a 5 permiten una 
inspección y comparación puntual o detallada de sus valores anuales, 
observándose una similitud notable tanto en sus valores anuales como en los 
tipos de sequías meteorológicas que definen, independientemente del método de 
estimación de la ETP. Lo anterior será ratificado de manera numérica en la Tabla 
6 y Tabla 7; en la Figura 2 se pueden observar los resultados de la Tabla 3, 
relativos a la estación Villa de Arriaga con k = 12 meses. 


Tabla 6. Evaluación de contrastes a través del ECM y del ESM entre la ETP y el 
RDIst anuales de la fórmula de Penman-Monteith, así como sus respectivos 
valores estimados con los tres métodos empíricos citados en las tres estaciones 
climatológicas indicadas del estado de San Luis Potosí, México (entre 
paréntesis, los valores mínimos de cada contraste). 


Estación: Villa de Arriaga Estación: Río Verde Estación: Xilitla 
k =3 k=6 k=12 k =3 k=6 k=12 k =3 k=6 k=12 


Concepto 
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ECM de la 
ETP anual 
de 
Thornth- 
waite 


668.3 


346.4 


110.8 


ECM de la 
ETP anual 
de Turc 


166.8 


94.2 


(39.1) 


ECM de la 
ETP anual 
de 
Hargrea- 
ves- 
Samani 


(89.1) 


(24.5) 


103.1 


ESM de la 
ETP anual 
de 
Thornth- 
waite 


662.8 


344.8 


104.8 


ESM de la 
ETP anual 
de Turc 


154.7 


(7.8) 


ESM de la 
ETP anual 
de 
Hargrea- 
ves- 
Samani 


(74.9) 


(-15.5) 


-97.5 


ECM de los 
RDIst 
anuales de 
Thornthwa 
¡te 


0.165 


0.210 


0.136 


0.072 


0.100 


0.097 


0.086 


0.117 


0.112 


ECM de los 
RDIst 
anuales de 
Turc 


0.062 


0.080 


0.062 


0.035 


0.044 


0.046 


0.082 


0.089 


0.108 


ECM de los 
RDIst 
anuales de 
Hargreave 
s-Samani 


(0.050) 


(0.062) 


(0.050) 


(0.029) 


(0.036) 


(0.039) 


(0.070) 


(0.076) 


(0.098) 


ESM de los 
RDlIst 
anuales de 
Thornth- 
waite 


3.104-1 
07 


ESM de los 
RDIst 
anuales de 
Turc 


ESM de los 
RDIst 
anuales de 
Hargrea- 
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ves- 
Samani 


Tabla 7. Severidad de las sequías meteorológicas, obtenida con el RDIst en las 
tres duraciones (k) estudiadas en meses, con cada uno de los criterios de 
estimación de la ETP aplicados en las tres estaciones climatológicas indicadas 
del estado de San Luis Potosí, México. 


Tipos Penman-Monteith Thornthwaite Turc Hargreaves-Samani 
a k =3 k= 6 k k =3 k=6 k k =3 k=6 k k =3 k=6 k= 
y 12 12 12 12 
meteor [N|[%|N|%|N|%|N %|N %|N|%|N|%|N|%|N|% | N|%|N|%|N|9% 
ológica | o lo) (o) (o) (o) (o) (o) (o) (o) (o) (o) (o) 
S 
(SMET 
dE 
Estación climatológica: Villa de Arriaga (NA/2 = 26.5) 
SMET 2|8|1|6/[1|5|2|8|1|6|1|6|1|8|1|6|1|5|1|8|1|6|1|5 
ligeras [| 0|3.|7|5.[|5|7.|2|4.|8|9.|8|4.|9|2.|7|3.|6|9.|9]|2.|7|3.|6|09. 
3 4 7 6 2 3 6 0 3 6 0 3 
SMET 1 2 1 2 2 2 2 2 
moder [2 |8.|5|9.|7|6.|2|7.|5|9.|6|1 2|8.|6|2.17|5.|2|8.|6]|2.|7]|5. 
adas 3 2 9 7 2 4 7 2 9 7 2 9 
SMET 1 1 1 1 1 1 1 
severa | 1|4.|3|1.|3|1 1/|3.|2|7.|3|0 1|4.|3/|1 3|(1.|1|4.|3|1.|3|1. 
S 2 5 5 8 7 7 3 1 1 3 1 1 
SMET 3. 3. 
extrem | 1/|4.[|1|3.|1|3.|1|3.|1/|/3.|1/|/3.[1|/4.[1/3.|1/7|1/|/4. 117 |11|3. 
as 2 8 8 8 8 6 3 7 3 7 
Estación climatológica: Río Verde (NA/2 = 27) 
SMET 1|6|2|7|1|6|1|6|1|7|1|6|1|6|2|7|1|6|1|6|2|7|1|6 
ligeras | 7 | 8.| 2 |8.|8|6.|8|9.|9|0.|9|7.|8|9.|2|8.|8|6.|8|9.|2|8.|8| 6. 
0 6 7 2 4 9 2 6 7 2 6 7 
SMET 1 1 2 2 2 2 2 L: 2 2 1 2 
moder [4|6.|4|4.[6|2.|6|3.|6|2.|6|1.|6|3.|4|4.|6|2.|6|3.|41|4.|6|2. 
adas 0 3 2 1 2 4 1 3 2 1 3 2 
SMET 7. 
severa | 2|8.|0|0.|2|7.|0|0.[|1|3.[|2|1|0|0.[1|3.|2|7.|0|0.|1|3.|2|7. 
s 0 0 4 0 7 0 6 4 0 6 4 
SMET 3. 
extrem |2|8.|2|7.|1|3.|2|7.|1|3.[1|6|2|7.|[1/3.|1/3.|2|7.|11/|3.|1|3. 
as 0 1 7 7 7 7 6 7 7 6 7 
Estación climatológica: Xilitla (NA/2 = 25) 
SMET 16|1|6|1|5|1|5|1|6|1|5|1|6|1|6J|1|5|1|6|1|5|1|5 
ligeras | 3|1.|4|0.|1|0.|2|7 5|0.|4|8.|4|3.|4/|0.|3|9.|3|1.|4|8.|3|9 
9 9 0 1 0 3 6 9 1 9 3 1 
SMET 1 1 3 1 2 2 1 1 1 1 2 1 
moder [3 |4.|4|7.|7|1.|4|9.|6|4.|6|5.|3|3.|4|7.|4/|8.1|3|4.|5|0.|41]8 
adas 3 4 8 0 0 0 6 4 2 3 8 2 
SMET 1 1 1 1 1 1 1 1 
severa [3 |4.|3|3.|2|9.|3|4 3|2.|2|8.|2|9.|3|3.|3|3 2|9.|3|2.|3|3 
S 3 0 1 3 0 3 1 0 6 5 5 6 
SMET 1 1 
extrem | 2|9.|2|8.|2|9.|2|9.|1|4.|2|8.|3|3 2|8.|2|9.|3|4.|2|8.|2|9 
as 5 7 1 5 0 3 6 7 1 3 3 1 
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(k = 12 meses) 


Figura 2. Comparación de los 53 valores del RDIst calculados con la ETP de 
Penman-Monteith (ordenadas) y de Hargreaves-Samani (abscisas), en la 
estación climatológica Villa de Arriaga, San Luis Potosí. 


Resultados del ECM y ESM 


En la Tabla 6 se han concentrado los valores numéricos del ECM y del ESM. El 
contraste anual de las estimaciones de la ETP indica que en los climas semiárido 
y templado-seco, el método de Thornthwaite es el que menos se aproxima a los 
resultados de la fórmula de Penman-Monteith, y el de Hargreaves-Samani es el 
más exacto. En el clima cálido-húmedo, los dos métodos anteriores se aproximan 
de manera casi igual, y el de Turc es el más exacto. El signo menos del ESM de 
los métodos de Turc y Hargreaves-Samani (última columna de la Tabla 6) indica 
que estos criterios sobreestiman la ETP con respecto a la de referencia. Los 
resultados anteriores definieron la manera en que se presentan los resultados de 
RDIs en la Tabla 3, Tabla 4 y Tabla 5. 
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El ECM relativo a los valores anuales del RDIst, en las tres duraciones analizadas, 
resultó mayor con el método de Thornthwaite y de orden de magnitud similar 
con los otros dos criterios, pero siempre condujo a un valor menor, en los tres 
climas estudiados, el método de Hargreaves-Samani. Con respecto a los valores 
del ESM obtenidos, en general son bajos, del mismo orden de magnitud y 
ratifican los resultados del ECM, variando su signo según el método de estimación 
de la ETP y la duración k. 


Severidad de las sequías meteorológicas (SMET) 


En la Tabla 7 se exponen las estimaciones relativas al número obtenido de cada 
tipo de SMET en cada una de las tres duraciones (k) y estaciones climatológicas 
procesadas. De manera general, para la duración de tres meses es donde se 
tienen las mayores dispersiones en el número de SMET, que en teoría debe ser 
igual a la mitad del número de años de registro (NA), valor que se indica para 
cada estación climatológica. Los porcentajes citados en la Tabla 7 se calcularon 
con el número encontrado de SMET, por ello suman 100%. 


Se observa que con el método de Thornthwaite se encuentran los porcentajes de 
cada tipo de SMET, que son más disímiles de los obtenidos con la ETP de 
referencia. Esto ocurre en las tres estaciones climatológicas, pero es más 
acusado en Río Verde. 


Se puede indicar que los porcentajes de cada SMET que definen los métodos de 
Turc y de Hargreaves-Samani son bastante similares a los obtenidos con la 
fórmula de Penman-Monteith. Lo anterior ratifica los resultados de la Tabla 6. 


Los valores numéricos de la Tabla 6 y Tabla 7 permiten concluir que no se tiene 
influencia significativa en los valores anuales de RDI« ni en los porcentajes de 
cada tipo de SMET que ellos definen, cuando se aplica el método de Hargreaves- 
Samani, en cualquiera de las tres estaciones climatológicas procesadas. El 
método de Thornthwaite resulta aplicable exclusivamente en el clima cálido- 
húmedo de la estación climatológica Xilitla. 


Conclusiones 


Los resultados de la aplicación del RDIst en las tres estaciones climatológicas 
procesadas del estado de San Luis Potosí, pertenecientes a climas diferentes, 
indican que no existe influencia significativa en los valores anuales del RDIst o 
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en sus porcentajes de cada tipo de sequía meteorológica que detectan cuando 
se emplean los métodos empíricos de Hargreaves-Samani y de Turc en la 
estimación de evapotranspiración potencial (ETP) mensual, en comparación con 
los resultados de la fórmula de Penman-Monteith tomada como referencia. 


Lo anterior permite establecer al RDIs como un índice de sequías meteorológicas 
robusto, prácticamente no dependiente del método de estimación de la ETP. 


El cálculo numérico de la ETP, según la fórmula de Penman-Monteith y el método 
de Hargreaves-Samani, tiene una diferencia notable en complejidad, y por ello 
el resultado de la Tabla 6, la cual indica que el ECM es el menor en los tres climas 
estudiados, tal criterio empírico es sumamente importante por su trascendencia 
práctica. 


Apéndice 1: fórmula de Penman-Monteith 


Ecuaciones teórica y operativa 


Penman, en 1948, fue el primero en obtener una ecuación que combina la 
energía requerida para sostener la evaporación y una descripción empírica del 
mecanismo de difusión por el cual la energía es removida de la superficie de 
evaporación, como vapor de agua (Shuttleworth, 1993). La fórmula de Penman 
originó un nuevo criterio de estimación de evaporación, denominado Método de 
Combinación. Varios investigadores modificaron la fórmula de Penman para 
tomar en cuenta los efectos de la evolución de las condiciones aerodinámicas con 
el crecimiento del cultivo, lo anterior a través de los factores de resistencia. Se 
distingue la resistencia de superficie (rs) que se ofrece al flujo de vapor de agua, 
en los estomas de las hojas, en la superficie del suelo y la resistencia 
aerodinámica (ra) que se debe a la fricción del flujo de aire sobre la superficie 
vegetal. Aunque los procesos de intercambio en la capa de vegetación son mucho 
más complejos, las mediciones y cálculos del flujo de calor latente AET han 
mostrado correlación alta, al menos para una superficie de pasto uniforme. Con 
tales modificaciones se llegó a la fórmula de Penman-Monteith teórica (Allen et 
al., 1998): 


_ A(Rn-G)+Pa' Cp (es-e)/Ta 
AET = A+y(1+T5/Ta) (A.1) 
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en la cual, AET es la velocidad de evapotranspiración en megajoule por m? y por 
día (MJ/m2/d); A, la pendiente en un punto de la curva de vapor de saturación 
contra la temperatura en kilopascal por *C (kPa/*C); Rn, la radiación solar neta 
en MJ/m?/d; G, el flujo de calor desde el suelo en MJ/m?/d; pa, la densidad media 
del aire a presión constante en kg/m; cp, el calor específico del aire a presión 
constante en MJ/kg/*C; (es - e), el déficit de presión de vapor del aire en kPa; 
Fa, la resistencia aerodinámica en s/m; y, la constante psicrométrica en kPa/*C, 
y rs es la resistencia de superficie en s/m. 


Al considerar una superficie hipotética de cultivo con 12 cm de altura, con una 
resistencia de superficie fija de 70 s/m, un albedo de 0.23, en crecimiento activo, 
que sombrea completamente al suelo y que no le falta agua, se obtiene la 
siguiente fórmula operativa de Penman-Monteith (Allen et a/., 1998): 


_0.408-4-(Rn—G)+y-[900/(Tt+273)]-uz"(es—e) (A.2) 


ET, A+ y:(1+0.34:uz) 


en donde ETo es la evapotranspiración de referencia en milímetros por día 
(mm/d) y los dos nuevos términos son los siguientes: Tt es la temperatura media 
del aire a dos metros de altura en *C y u» es la velocidad media del viento a dos 
metros de altura en m/s. El valor de rs = 70 s/m corresponde a una superficie 
del suelo moderadamente seca, resultante de una frecuencia semanal de riego. 


Para el paso de la ecuación (A.1) a la (A.2), se toma en cuenta que una lámina 
de agua en mm/d se puede expresar en términos de la energía recibida por 
unidad de área. Tal energía se refiere al calor necesario para evaporar tal espesor 
de agua y se conoce como calor latente de evaporación (A), es una función de la 
temperatura del agua (Ta) y se calcula con la siguiente ecuación en MJ/kg (Allen 
et al., 1998): 


2 = 2.501 — 0.002361 «Ta (A.3) 


Como el valor de A no cambia mucho con Ta, se utiliza Ta = 20 *C y entonces A 
es cerca de 2.45 MJ/kg, es decir, se requieren 2.45 MJ para evaporar un 
kilogramo de agua o un litro; entonces, una lámina de un mm de agua es 
equivalente a 2.45 MJ/m?, ya que 1 mm por m? es un decímetro cúbico, es decir, 
un litro. El primer coeficiente numérico de la ecuación (A.2) convierte la radiación 
expresada en MJ/m?/d a evaporación en mm/d y equivale al valor inverso de A 
(1/4 = 0.408). 


La ecuación (A.2) se puede aplicar en intervalos de un día, 10 días, un mes o 
incluso la duración total de crecimiento del cultivo o un año. Para obtener ETo 
en mm/h se cambia el numerador del paréntesis rectangular a 37 y todas las 
variables son por hora no por día. Con propósitos de verificación de sus 
resultados, se indica que en las regiones tropicales húmedas con temperatura 
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media (Tt) moderada (=20 *C) ETo varía de 3 a 5 mm/d y con Tt caliente (> 
30%C) oscila de 5 a 7 mm/d; aumentado estos intervalos una unidad en las zonas 
áridas (Allen et a/., 1998). En México, ya se han realizado aplicaciones de la 
ecuación (A.2) por González-Camacho, Cervantes-Osornio, Ojeda-Bustamante, 
y López-Cruz (2008), y Chávez-Ramírez et a/. (2013). 


Estimación de los parámetros A y y 


Todas las expresiones que se presentan a continuación proceden de Allen et al. 
(1998) y se orientan a estimar la evapotranspiración potencial (ETP'¡) del mes j¡ 
de cada año /, que es requerida en la aplicación de la ecuación (1). La pendiente 
(A en kPa/*C) en la curva de presión de vapor de saturación en un punto relativo 
a la temperatura media del aire (Tt) en *C se calcula con la expresión: 


17.27.18; 
4 089:|0.6108:exp | ——— 
Ttl4237.3 


Ay S a di 


La constante psicrométrica (y en kPa/*“C) se determina con la expresión 
siguiente: 


i_ CpP _ 1.62862:10"3-P A 
Y EA; x; (043) 


en donde cp = 1.013-103 MJ/kg/*C y e = 0.622 es el cociente del peso molecular 
del vapor de agua al del aire; A se estima con la ecuación (A.3) para el valor de 
Tt¿ en “C, y P es la presión atmosférica del lugar en kPa; se estima con la 


ecuación: 


¿526 
P=101.3: (E) (A.6) 
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donde z es la altitud en metros sobre el nivel del mar. 


Estimación de las radiaciones Rn y G 
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La radiación neta (Rn) equivale a la diferencia entre la radiación solar neta de 
onda corta incidente (Rns) y la radiación solar neta de onda larga que se emite 
o libera (Rn)), es decir: 


Rn = Rns — Rnl (A.7) 


La Rns es la diferencia entre la radiación solar incidente (Rs) y la reflejada, por 
lo cual se estima con la expresión: 


Rnsj = (1 — a)-Rsj(A.8) 


en donde a es el albedo o coeficiente de reflexión de la cubierta vegetal, 
adimensional. Un valor de 0.23 se adopta para el cultivo hipotético del pasto de 
referencia. Rs; se debe expresar en MJ/m2/d, por lo cual los valores promedio 
mensual en cal/cm?/d leídos en los mapas propuestos por Almanza y López 
(1975), o bien por Hernández et a/. (1991) se deben multiplicar por 0.041868 
para obtener MJ/mY/d. 


La velocidad de emisión de energía de onda larga es proporcional a la 
temperatura absoluta de la superficie elevada a la cuarta potencia; relación que 
se conoce como Ley de Stefan-Boltzmann. Como el vapor de agua, las nubes, el 
dióxido de carbono y el polvo absorben y emiten radiación de onda larga, su 
balance o flujo neto que deja la superficie terrestre se estima corrigiendo la ley 
de Stefan-Boltzmann por la humedad relativa y la nubosidad, según la ecuación 
siguiente: 


Rnii o(rep)' : (0.3 E 0.14. [e)) , [Las ; (5) me 0.35| (A.9) 


en la cual, O = 4.903-10"? MJ/K*/m?/d es la constante de Stefan-Boltzmann; Tt;, 
la temperatura media del mes en grados Kelvin, igual a los grados centígrados 
(9C) más 273.16; ej la presión parcial de vapor actual en kPa, y Rso es la 
radiación solar en días despejados o sin nubosidad en MJ/m2?/d; se estima con la 
expresión: 


Rso = (0.75 + 2-10—5-z)-Rej(A.10) 
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en donde z nuevamente es la altitud del lugar y Re; es la llamada radiación 
extraterrestre en MJ/m?/d. El cociente Rs/Rso debe ser menor que la unidad. Las 
estimaciones de ej de Rej y de los dos términos faltantes (es y u2) de la ecuación 
(A.2) se detallan en el Apéndice 2. 


Al tomar en cuenta que el flujo de calor desde el suelo (G) es reducido en 
comparación con Rn, se emplea una expresión muy simple para su estimación, 
la cual acepta que la temperatura del suelo es semejante a la del aire: 


TEj+TEj- 2 


G = Sa 


Az (A.11) 


donde cs = 2.10 MJ/m?/*C es la capacidad calórica del suelo; Ad es el intervalo 
en días, y Az es la profundidad afectada del suelo; para lapsos de un mes o 
mayores se considera igual a dos metros. Con base en tales valores numéricos, 
la ecuación (A.11) del primer mes, subsiguientes y último son las siguientes: 


Gj=1 = 0.14 - (Tt;4, — Tt;) (A.12) 


Gj=12NA = 0.14 , (Tt; e Tt;_1) (A.14) 


En la ecuación (A.14), NA es el número años del registro climático procesado. 


Apéndice 2: estimaciones climáticas 
complementarias 


Radiación extraterrestre 


Rey es la radiación solar en el tope de la atmósfera en cal/cm?/d; está tabulada 
de forma mensual, y es función de la latitud del lugar (q) en grados. Para evitar 
la interpolación de Rey se desarrollaron 12 polinomios de Newton de tercer grado, 
cuya fórmula es aplicable a latitudes de 10 a 40 grados norte (Campos-Aranda, 
2005a): 
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Re; = bo + b(p — 10) + b2(p — 10)(p — 20) + ba(p — 10)(p — 20)(p — 30) (A.15) 


Los coeficientes b; se tienen en la tabulación siguiente: 


Mes | bo | bi b2 b3 Mes | bo bi b2 b3 

Ene |760 -12 -0.075| 1/600 | JUL | 880 5 |-0.100 | -1/1200 
Feb [820 | -9 |-0,.100 | 1/1200 | AGO | 890 2 |-0.125 | -1/1200 
Mar |875| -5 |-0.125| 1/1200 | SEP | 880 | -2.5 | -0.150 | 1/1200 
Abr [895 | O |-0.125 | -1/1200 | OCT |840| -8 |-0.075 | -1/1200 
May |890| 4 |-0.100 | -1/400 NOV | 780 | -11.5 | -0.025 | -1/300 
Jun |875| 6 |-0.100 | -1/600 | DIC | 740 | -12.5 | -0.075 | 1/1200 


Los valores estimados de Re, con la ecuación (A.15) se deben multiplicar por 
0.041868 para obtenerlos en MJ/m?/d. 


Presión parcial de vapor 


Para estimar e'; de la expresión (A.9), se debe recordar que el aire entre más 
caliente puede contener mayor cantidad de vapor de agua, cuyo máximo es la 
presión parcial de vapor de saturación (es) en kilo Pascales (kPa), la cual es 
función de la temperatura y se estima con la expresión siguiente: 


17.27-Tt; 


esj = 0.6108- exp | (A.16) 


(Ttj+237.3) 


Cuando existe menos cantidad de vapor de agua que el máximo, tal presión 
parcial de vapor de agua se designa por e y entonces la humedad relativa (HR) 
en porcentaje es: 


HR = o 100 (A.17) 


Si el aire conteniendo una cantidad de vapor de agua igual a e se va enfriando, 
se llega a un punto en que e se convierte en es y esa temperatura se denomina 
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de punto de rocío (t«*), misma que se estima con la ecuación (A.16), 
despejándola; la cual conduce a: 


2 237.3-[m(et/0.6108)] 
yo 17.27-[Im(e¿/0.6108)| 


(A.18) 


Habiendo observado que t* es muy próxima a la temperatura promedio mensual 
de mínimas (t';), se ha establecido una manera simple de estimar el valor de e';. 
Lo anterior fue verificado hace casi tres décadas por Arteaga-Ramírez (1989), y 
recientemente por Cervantes-Osornio, Arteaga-Ramírez, Vázquez-Peña, Ojeda- 
Bustamante y Quevedo-Nolasco (2013). 


En la Tabla Al se han concentrado los datos promedio mensuales de humedad 
relativa (ecuación (A.17)), presión parcial de vapor (e) y temperatura mínima 
(t) de cinco observatorios meteorológicos (SARH, 1982), circundantes a las tres 
estaciones climatológicas que serán procesadas. Con base en la ecuación (A.18), 
usando 6.108 en lugar de 0.6108, ya que e está en mbar, se obtuvieron las 
temperaturas de punto de rocío (t*) correspondientes. Después se calcularon las 
diferencias entre t y t* para realizar su inspección y establecer las correcciones 
al valor de t, pues en teoría tales diferencias deberían ser cercanas a cero. 


Tabla A1. Valores promedio mensuales de varios elementos climáticos en los 
cinco observatorios meteorológicos indicados. *Precipitación media anual. 


Descripción Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic Anual 
Observatorio meteorológico: Saltillo (Coah.); PMA* = 269.4 mm 

Humedad relativa 62 59 54 54 58 62 65 68 72 70 64 62 62 

(%) 

Tensión de vapor 9.1 |96| 90 | 13.0 | 15.6 | 17. 17. 17.5 | 16.6 | 13. 11.5 | 9.1 13.3 

(mbar) 9 6 9 

Punto de rocío (t* 2C) | 5.6 | 6.4 5.4 | 10.9 | 13.6 | 15. 15. 15.4 | 14.6 | 11. 9.0 5.6 11.2 
8 5 9 

Temp. mínima (t *C) 5.2 |6.8 | 8.7 | 12.7 | 14.7 | 16. 16. 16.2 | 14.5 | 11. 8.0 6.3 11.4 
4 5 6 

Diferencias t - t«+ (9C) S 0.4 3.3 1.8 1.1 0.6 1.0 0.8 -0,1 | -0.3 | -1.0 | 0.7 0.2 

0.4 
Observatorio meteorológico: San Luis Potosí (SLP). PMA = 315.4 mm 

Humedad relativa 51 43 39 37 47 56 60 61 65 63 57 56 52 

(%) 

Tensión de vapor 75 |6.8 | 7.3 7.9 11.3 | 13. 13. 13.3 | 13.8 | 12. 9.6 8.3 10.4 

(mbar) 4 8 0 

Punto de rocío (t* *C) | 2.9 | 1.5 2.5 3.6 8.8 11. 11. 11.2 | 11.8 | 9.7 6.4 4,3 7.5 
3 8 

Temp. mínima (t *C) 6.2 | 74 | 99 | 119 | 134 | 14. 13. 13.5 | 13.2 | 10. 8.1 6.5 10.7 
2 5 7 
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Diferencias t - t* (9C) | 3.3 | 5.9 | 7.4 8.3 4.6 2.9 1.7 2.3 1.4 1.0 1.7 2.2 3.2 
Observatorio meteorológico: Río Verde (SLP). PMA = 484.9 mm 
Humedad relativa 73 70 64 64 66 69 73 72 77 76 76 75 71 
(%) 
Tensión de vapor 12. | 13. | 148 | 17.2 | 19.7 | 21. 20. 21.4 | 20.9 | 18. | 15.6 | 13. 17.5 
(mbar) 8 6 3 9 4 7 
Punto de rocío (t* *C) | 10. | 11. | 12.8 | 15.1 | 17.3 | 18. 18. 18.6 | 18.2 | 16. | 13.6 | 11. 15.4 
6 5 5 2 2 6 
Temp. mínima (t *C) 9.2 | 10. | 12.9 | 16.0 | 17.9 | 19. 18. 18.3 | 17.3 | 15. | 122 | 9.6 | 14.7 
8 0 2 0 
Diferencias t - t* (9C) 2 > 0.1 0.9 0.6 0.5 0.0 -0.3 -0.9 | -1.2 | -1.4 - -0.7 
1.4 | 0.7 2.0 
Observatorio meteorológico: Aguascalientes (Ags.). PMA = 537.2 mm 
Humedad relativa 57 52 46 43 46 59 65 67 69 64 59 61 57 
(%) 
Tensión de vapor 8.8 |8.9| 9.4 | 10.1 | 12.2 | 14. 15. 15.6 | 15.2 | 13. | 10.5 | 9.6 | 11.9 
(mbar) 8 3 0 
Punto de rocío (t* *C) | 5.1 | 5.3 | 6.1 7.1 9.9 12. 13. 13.6 | 13.2 | 10. 7.7 6.4 9.5 
8 3 9 
Temp. mínima (t *C) 4.6 | 6.1 8.2 | 11.0 | 13.4 | 14. 14. 13.9 | 13.3 | 10. 7.2 5.5 10.2 
8 1 6 
Diferencias t - tx (9C) 5 0.8 | 2.1 3.9 3.5 2.0 0.8 0.3 0.1 | -0.3 | -0.5 - 0.7 
0.5 0.9 
Observatorio meteorológico: Tampico (Tam.). PMA = 985.9 mm 
Humedad relativa 81 81 80 82 81 82 80 80 81 79 79 80 80 
(%) 
Tensión de vapor 17. | 19. | 21.1 | 25.5 | 28.5 | 30. 30. 30.4 | 29.7 | 25. | 21.5 | 18. 24.9 
(mbar) 7 3 5 3 9 9 
Punto de rocío (t* *C) | 15. | 16. | 18.4 | 21.4 | 23.2 | 24. 24. 24.3 | 23.9 | 21. | 18.7 | 16. 21.0 
6 9 4 3 7 6 
Temp. mínima (t *C) 14. | 15. | 17.6 | 20.8 | 22.9 | 23. 23. 24.1 | 23.1 | 21. | 18.3 | 15. 20.1 
1 7 9 8 2 7 
Diferencias t - t«* (9C) a 5 -0.8 | -0.6 | -0.3 | -0.5 | -0.5 | -0.2 -0.8 | -0.5 | -0.4 S -0.9 
1.5 | 1.2 0.9 


Para la estación Villa de Arriaga se combinan las correcciones en *C sugeridas 
en Saltillo (de septiembre a febrero con cero, y de marzo a agosto con -1.00) y 
Aguascalientes (de julio a febrero con cero, y de marzo a junio con -2.00), 
proponiéndose lo siguiente: septiembre a febrero con cero y marzo a agosto con 
-1.50. No se utilizan los datos del observatorio de San Luis Potosí, por no 
considerarse confiables. En el observatorio de Río Verde, las correcciones 
recomendadas son de marzo a agosto cero y de septiembre a febrero con +1.20. 
Por último, en la estación Xilitla se aplicará la corrección sugerida para el 
observatorio de Tampico, que es sumar en todos los meses 0.60 *C. 


Velocidad media del viento 
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Por último, con respecto a la estimación de la velocidad media del viento (v) a 
utilizar en la ecuación (A.2), la FAO ha establecido cuatro valores según el tipo 
de vientos que ocurren en la región: (1) débil, de 0.50 a 1.0 m/s; (2) moderado, 
de 1 a 3 m/s; (3) severo, de 3 a 5 m/s, y (4) fuerte > 5 m/s. También señala 
que la velocidad promedio a 2 m de altura en más de 2 000 estaciones 
meteorológicas del mundo es de 2 m/s. 


Campos-Aranda (2005b) procesó datos de velocidad media del viento (m/s) 
mensual en 31 observatorios meteorológicos que cuentan con tales datos, en 
amplitudes de registro que variaron de 8 a 20 años. Encuentra que la moda de 
2.05 m/s verifica el valor medio mundial observado y obtiene una mediana 
muestral y poblacional del orden de 2.3 m/s; recomendado tal valor como 
promedio mensual. 


Radiación solar con datos de temperatura 


Cuando no se disponga de los mapas de Rs; de Almanza y López (1975) o de 
Hernández et al. (1991), o bien no se quiera emplear valores promedio 
mensuales, se puede realizar una estimación con base en la diferencia mensual 
entre las temperaturas máxima (7) y mínima (t) del aire, ya que tal diferencia 
está relacionada con el grado de nubosidad de la localidad. En general, en los 
días despejados o sin nubosidad se generan altas temperaturas durante el día y 
bajas por la noche debido a que no se absorbe ni se regresa la radiación de onda 
larga. Lo contrario ocurre en días nubosos. La ecuación empírica operativa es 
(Allen et al., 1998): 


Rs; = ka: Re; [T/ —tj (A.19) 


en la cual Rs; es la radiación solar en MJ/m?/d; ka, un factor de ajuste que varía 
de 0.16 a 0.19 con unidades 1/*C, y Re, es la radiación extraterrestre en 
MJ/m2/d. El valor de ka = 0.16 se emplea en localidades del interior, donde las 
masas de aire no son influenciadas por el mar y el de ka = 0.19 se utiliza en 
zonas costeras donde sí hay tal afectación. Allen et al. (1998) también indican 
cómo se pueden transportar a la localidad en estudio datos de Rsj de una estación 
meteorológica cercana. 
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Apéndice 3: criterios empíricos de estimación de la 
ETP 


Método de Thornthwaite 


Estima la evapotranspiración potencial (ETP'¡) del mes j del año ¡ en función 
exclusivamente de la temperatura media (Tt'¡) en *C y de la latitud del lugar 
(LAT) en grados, su ecuación es (Mather, 1977; Campos-Aranda, 2005a; Xu, 
Singh, Chen, 8 Chen, 2008): 


10-Tél 


ETP¡ =16: ( a) "Fc; (A.20) 


en la cual IC; es un índice de calor anual, igual a la suma de los 12 mensuales, 
que son: 


El exponente a es función del IC;, con la siguiente ecuación empírica: 
ar = 0.4924 + 1.792 : 107? + IC; — 7.71: 107? -IC2 + 6.75: 107 - IC? (A.22) 


Por último, Fc; es un factor correctivo promedio mensual función de la latitud del 
lugar y del número de días del mes (ndm); su fórmula es: 


Fo =(2)-(53) (4.23) 


30 


donde N es el soleamiento máximo o número máximo de horas con sol promedio 
mensuales. Para su estimación en la república mexicana, Campos-Aranda 
(2005a) propuso la siguiente expresión empírica: 


N = A + B [sen (30 nm + 83.5)] (A.24) 
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donde nm es el número de mes, con 1 para enero y 12 para diciembre; A y B 
son constantes función de la latitud del lugar (LAT) en grados, con las 
expresiones siguientes: 


> 
ll 


12.09086 + 0.00266 - LAT (A.25) 


B = 0.2194 — 0.06988 - LAT (A.26) 


Método de Turc 


A inicios de la década de 1960, Turc propuso la ecuación siguiente para estimar 
la ETP'; mensual y de 10 días, la cual es función de la temperatura media mensual 
Tt!,, de la radiación solar incidente Rs; expresada en cal/cm?/d, y de la humedad 
relativa media mensual (Turc, 1961; Campos-Aranda, 2005a; Xu et al., 2008): 


ETP) =cj ES +50): FC (A.27) 


e Tt; ze 


El coeficiente c toma los valores siguientes: 0.40 para meses de 30 o 31 días, 
0.37 para febrero y 0.13, para un lapso de 10 días. El factor correctivo se aplica 
cuando la humedad relativa media mensual (HR!) es inferior a 50%; su expresión 
es: 


FCi=1 + 11) (4.28) 


Método de Hargreaves-Samani 


Propuesto a inicios de la década de 1980, estima la evapotranspiración potencial 
(ETP'¡) media diaria en milímetros, sólo a partir de la temperatura media (Tt';) 
expresada en grados Fahrenheit y de la radiación solar incidente media diaria 
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(Rs'j) indicada en milímetros de lámina de agua evaporada; su ecuación es 
(Hargreaves € Samani, 1982; Campos-Aranda, 2005a; Xu et al., 2008): 


ETP¡ = 0.0075 - Rs! -Tt! (A.29) 


La radiación solar incidente (Rs';) se puede estimar con la fórmula de Angstróm 
(Jáuregui-Ostos, 1978; Allen et al., 1998) cuando se tienen datos de insolación 
o soleamiento real (n), o bien con los mapas mensuales disponibles para la 
república mexicana (Almanza 8 López, 1975; Hernández et al., 1991), que la 
reportan en cal/cm?/d. Para la transformación de Rs; a lámina de agua evaporada 
por día se emplea la fórmula siguiente: 


Rsi = + (A.30) 
J 


J H 


donde Hv'; es el llamado calor latente de evaporación o energía necesaria en 
calorías para evaporar 1 g o un cm3 de agua; se estima con la expresión 
siguiente, estando la temperatura media (T+;) mensual en *C: 


Hvj = 595.9 — 0.55: Tt¡ (A.31) 


En Campos-Aranda (2005a) se puede consultar la otra fórmula empírica de 
Hargreaves-Samani, que es función de la radiación solar extraterrestre, y de las 
temperaturas media y mínima mensuales. Una aplicación de tal criterio se 
expone en Campos-Aranda (2014). 
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